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We investigate the assembly of groups and clusters of galaxies using the Millennium 
,i_ > , dark matter simulation and the associated Millennium gas simulations, and semi-analytic 

catalogues of galaxies. In particular, in order to find an observable quantity that could be 
used to identify early-formed groups, we study the development of the difference in magni- 
tude between their brightest galaxies to assess the use of magnitude gaps as possible indi- 
cators. We select galaxy groups and clusters at redshift z = 1 with dark matter halo mass 
-W(i?20o) ^ 10 13 /i -1 Mq, and trace their properties until the present time (z — 0). Fur- 
ther constraints are applied to keep those galaxy systems for which the X-ray luminosity 
^x.boi ^ 0.25 x 10 42 h~ 2 erg s~ x At redshift z = 0. While it is true that a large magnitude gap 
between the two brightest galaxies of a particular group often indicates that a large fraction of 
its mass was assembled at an early epoch, it is not a necessary condition. More than 90% of 
fossil groups defined on the basis of their magnitude gaps (at any epoch between < z < 1 ) 
cease to be fossils within 4 Gyr, mostly because other massive galaxies are assembled within 
their cores, even though most of the mass in their haloes might have been assembled at early 
times. We show that compared to the conventional definition of fossil galaxy groups based 
on the magnitude gap Ami2 ^ 2 (in the i?-band, within 0.5 i?200 of the centre of the group), 
an alternative criterion A77114 ^ 2.5 (within the same radius) finds 50% more early-formed 
systems, and those that on average retain their fossil phase longer. However, the conventional 
■ criterion performs marginally better at finding early-formed groups at the high-mass end of 

groups. Nevertheless, both criteria fail to identify a majority of the early-formed systems. 

Key words: cosmology: theory — galaxies: formation — galaxies: kinematics and dynamics 
— hydrodynamics — methods: numerical 



1 INTRODUCTION 

Although existing observations of the large scale structure of the 
Universe overwhelmingly favour cold dark matter cosmologies 
with hierarchical structure formation, the paradigm faces chal- 
lenges both from the existence of luminous passive galaxies at high 
redshift, and the abundance of low-mass g alaxies in the local uni- 
verse (e.g Baugh 2006; Bal ogh et alj|2008h . Galaxies dominate the 
visible universe and any cosmological model is expected to repro- 
duce the observed global properties of galaxies, at least statistically, 
in the first instance. A significant fraction of the evolutionary life 
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of many galaxies is spent in the environment of small systems (i.e. 
groups), where close interactions and mergers of galaxies occur 
with hig her efficiency than in massive haloes such as galaxy clus- 
ters (e.g. lMiles et alj|2004h . The observable properties of the bary- 
onic content of a group, which consists of the constituent galaxies 
and the inter-galactic medium (IGM), should be linked to mass as- 
sembly of the host group and its subsequent evolution. Among such 
observable properties, the "magnitude gap", i.e., the difference in 
the magnitudes of the two brightest galaxies, has been widely used 
as an optical parameter related to mass assembly of groups and 
clusters. 

Several studies show that a system of galaxies, where 
most of the mass has been assembled very early, de- 
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velops a larger magnitude gap compared to systems that 
form later. The idea is supported in observational samples 
(e.g. [P onman et al. 1994; Khosrosha hi. Jones. & Ponmanl |2004 
Khosroshah i. Ponman. & Jonesl |2007|). theoretical studies (e.g. 
Milosavrievic et alj 120061 ; Ivan den Bosch et all |2007|) and in de- 
tailed analysis of N-body numerical sim ulations " barnesl [l989l : 
iD'Onghia et alj|2005l ; lbariush et ai]|2007l) . These studies also pre- 
dict that such early-formed galaxy groups or clusters should be re- 
laxed and relatively more isolated systems in comparison to their 
later- formed counterpar ts. 

Ijones et alj (2003) defined such early-formed systems (also 
known as fossils) to have a minimum X-ray luminosity of Lx.bol ^ 
0.25 x 10 42 /i _2 erg s _1 , and a large magnitude gap in the i?-band 
between their first two brightest galaxies, i.e. Ami2 ^ 2.0, to 
distinguish them from late-formed groups and clusters. Recent re- 
search based on the Sloan Digita l Sky Survey has either used only 
the optical criterion dSantos. Mendes de Oliveira, & Sodre 120071) 
or both optical and X - ray criteri a (Eigent haler & Zeilingedl2009l ; 
IVoevodkin et al J 120091 ; lLaBarberaetal.ll2009h to identify bssils. 
In the former case, if the optical definition is solely employed, the 
chance of identifying truly early-formed systems diminishes, since 
a large fraction of the systems detected might be in the stage of 
collapsing for the first time, and so would not be X-ray luminous. 

From numerical simulations, Ivon Benda-B eckmann et all 

J2008I) found that Ami2 may not be a good indicator for iden- 
tifying early-formed groups, since the condition would no longer 
be fulfilled when a galaxy of intermediate magnitude fell into the 
group. This study was based on simulations of dark matter particles 
only, and does reveal how frequently such a situation would arise. 
Other, potentially more robust, magnitude gap criter i a have been 
considered in the literature. For example. ISales et al] (2007) finds 
that the difference in magnitude between the first and 10 th brightest 
galaxies in three fossil groups span a range of ~3-5, in agreement 
with their results from the analysis of the Millennium d ata together 
with the semi-analytic catalogue of ICroton etai] J2006l) . 

The overall number of observed fossil galaxy groups 
is small, making it difficult for observed systems to be sta- 
tistically compared to simulated systems. In spite of their 
low space density, fossils have been used to tes t mod- 
els of cosmological evolution (iMilosavlievic et al] 120061: 
Khosroshahi, Ponman, & Jones 2007; von Benda-Beckmann et al. 
20081 ; iDfaz-Gimenez. Muriel. & Mendes de Oliveiral 120081) . since 
the criteria for observationally identifying such systems are simple 
(Jones et al. 20031; [Khosrosh ahi. Ponman, & Jones 200(1), and it is 
generally assumed that they are the archetypal relaxed systems, 
consisting of a group-scale X-ray halo, the optica l image being 
domi nated by a giant elliptical galaxy at the core JPonman et al] 
Il994h . Indeed, if they are relaxed early-formed systems, fossil 
groups can be the ideal systems in which to study mechanisms 
of feedback, and the interaction of central AGN and the IGM of 
the group, since the effect of the AGN would not be compl icated 
by the effect of recent mergers djetha et al] 120081 2009). The 
lack of recent merging activity would also predict the absence of 
current or recent star for mation in early-type galaxies be l ongin g 
to fossil groups (e.g. iNolan. Ravchaudhurv. & Kabanl 120071) . 
and the relative dearth of red star-forming galaxies compared 
to similar elliptical-dominated non-fossil groups and clusters 
dMahaian & Ravchaudhurvj |2009 ) . 

The use of cosmological simulations in the study of the evolu- 
tion of galaxy groups will have to employ a semi-analytic scheme 
for simulating galaxies, and the results will be dependent on the 
appropriate characterisation of the models that describe galaxy for- 



mation and evolution. Once the hierarchical buildup of dark mat- 
ter haloes is computed from N-body simulations, galaxy formation 
is modelled by considering the rate at which gas can cool within 
these haloes. This involves assumptions for the rate of galaxy merg- 
ing (driven by dynamical friction) and the rate and efficiency of 
star formation and the associated feed back in individual galaxies 
dCroton et alj|2006l ; lBower et al]|2006b 

In a previous paper dPariush et ai] |2007). we studied the for- 
mation of fossil groups in the Millennium Simulations, and showed 
that the conventional definition of fossils (namely a large magni- 
tude gap between the two brightest galaxies within half a virial ra- 
dius and a lower limit to the X-ray luminosity Lx.boi ^ 2.5 x 
10 42 ft~ 2 erg s _1 ), results in the identification of haloes which are 
~ 10%-20% more massive than the rest of the population of 
galaxy groups with the same halo mass and X-ray luminosity, when 
the Universe was half of its current age. This clearly indicates an 
early formation epoch for fossils. In addition, it was shown that 
the conventional fossil selection criteria filter out spurious systems, 
and therefore there is a very small probability for a large magni- 
tude gap in a halo to occur at random. The fraction of late-formed 
systems that are spuriously identified as fos sils was found to be 
~ 4-8%, almos t independent of halo mass JDariush et al]|2007l; 
ISmith et al.l2010l) . Another important outcome of the lDariush et all 
d2007t) study was the consistency between the space density of fos- 
sils found in the simulations and that from observational samples. 

Although the results from this previous analysis of the Millen- 
nium simulations are shown to be in fair agreement with observa- 
tion, we did not investigate the evolution of the magnitude gap in 
either fossil or control groups with redshift. Furthermore, the num- 
ber of (fossil or control) groups used was small (~400), which did 
not allow us to fully explore the connection between the halo mass 
and magnitude gap in such systems. 

In this paper, we select early-formed galaxy groups from the 
Millennium simulations, purely on the basis of their halo mass evo- 
lution from present time up to redshift z « 1.0, and, with the help 
of associated semi-analytic catalogues, study the evolution of the 
magnitude gap between their brightest galaxies. Our aim is to (a) 
investigate how well the conventional optical selection criterion, 
namely the Am ^ 2 gap between the two brightest galaxies, is able 
to identify early-formed galaxy groups, and (b) to find whether we 
can find a better criterion to identify groups that have assembled 
most of their mass at an early epoch. 

In §2, we describe the various simulation suites used in this 
work, and in §3 the data we extract from them. In §4, we study 
in detail the evolution with epoch of various measurable parame- 
ters for a large sample of early-formed "fossil systems" and two 
comparable sample of control systems, and compare these prop- 
erties. In §5, we examine the case for a revision of the criteria to 
observationally find fossil in order to ensure a higher incidence of 
genuine early-formed systems. We summarise our conclusions in 
§6. We adopt Ho = 100 ft km s 1 Mpc 1 for the Hubble constant, 
with h = 0.73. 



2 DESCRIPTION OF THE SIMULATIONS 
2.1 The Millennium Simulation 

The Millennium Run consists of a simulation, in a Universe con- 
sistent with concordance ACDM cosmology, of 2160 3 particles of 
individual mass 8.6 x 10 s h' 1 M , within a co-moving periodic 
box of side 500/i _1 Mpc, employing a gravitational softening of 



5h 1 kpc, from redshift z = 127 to the present day (Sprin ge! et al.l 
120051) . The basic setup is that of an inflationary Universe, domi- 
nated by dark matter particles, leading to a bottom-up hierarchy 
of structure formation, which involves the collapse and merger 
of small dense haloes at high redshifts, into the modern-day ob- 
served large virialised systems such as groups and clusters. The 
cosmological parameters used by the Millennium Simulation were 
fi A = 0.75, Q M = 0.25, D. b = 0.045, n = 1, and a 8 = 0.9, and 
the Hubble parameter h = 0.73. 

Dark matter haloes are found in this simulation down to a 
resolution limit of 20 particles, yielding a minimum halo mass of 
1.72x 10 10 /i _1 Mq. Haloes in the simulation are found using a 
friends-of-friends (FOF) group finder, configured to extract haloes 
with overdensities o f at least 200 relative to the critical density 
dSpringeletai]|2005h . Within a FOF halo, substructures or sub- 
haloes are ide n tified using the SUBFIND algorithm developed by 
Sprin geTet al.l ( 1200 lh . and the treatment of the orbital decay of 
satellites is described in the next section. 

During the Millennium Simulation, 64 time-slices of the loca- 
tions and velocities of all the particles were stored, spread approx- 
imately logarithmically in time between z = 127 and z = 0. From 
these time-slices, merger trees were built by combining the tables 
of all haloes found at any given output epoch, thus enabling us to 
trace the growth of haloes and their subhaloes through time within 
the simulation. 



2.2 Semi-analytic galaxy catalogues 

2.2.7 The Croton et al. semi-analytic catalogue 

Croton ct al. (2006) simulated the growth of galaxies, and their 
central supermassive black holes, by self-consistently implement- 
ing sejruj^maly^ic mode ls of galaxies on the dark matter haloes of 
fhe lSpringel et al.U2005h simulation, Their semi-analytic catalogue 
contains 9 million galaxies at z — brighter than absolute mag- 
nitude Mr — 5 log h = —16.6, "observed" in B, V, R, I and K 
filters. The models focus on the growth of black holes and AGN 
as sources of feedback sources. The inclusion of AGN feedback in 
the semi-analytic model (allowing central cooling to be suppressed 
in massive haloes that undergo quasi-static cooling), and its good 
agreement with the observed galaxy luminosity function, distribu- 
tion of galaxy colours and of the clustering properties of galaxies, 
make this catalogue suitable for our study. 

In this semi-analytic formulation, galaxies initially form 
within small dark matter haloes. Such a halo may fall into a larger 
halo as the simulation evolves. The "galaxy" within this halo then 
becomes a satellite galaxy within the main halo, and follows the 
track of its original dark matter halo (now a subhalo), until the 
mass of the subhalo drops below 1.72x 10 10 /i -1 Mq. This limit 
corresponds to the 20-particle limit for dark haloes in the original 
Millennium Simulation. At this point the galaxy is assumed to spi- 
ral into the centre of the halo, on some fraction of the dynamical 
friction tim escale, where it me rges with the central galaxy of the 
larger halo dCroton et alj2006l) . 



2.2.2 The Bower et al. semi-analytic catalogue 

The lBower et al. model also makes use of the Millennium 

Simulation, but utilises merger tr ees constructed with the algorithm 
described bv lHarker et al.l d2007l) . The cooling of gas and the sub- 
sequent formation of galaxies and black holes is followed through 



The mass assembly of galaxy groups 3 

the merging hier archy of each tree utilising the GALFORM semi- 
analy tic model dCole et alj|"2000l: iBower et al.ll2006l : lMaibon et al.l 
120071) . At z = this results in 4,491,139 galaxies brighter than a 

limiting absolute magnitude of Mr — 5 log h = —19.4. 

In addition to feedback from supernovae, the IBower et al.l 

d2006l) model accounts for energy input from AGN, resulting in 
a suppression of cooling in the hot atmospheres of massive haloes. 
The resulting galaxy population is in excellent agreement with the 
observed z = galaxy luminosity function in B and K bands, the 
z — colour distribution and also with the evolution of the galaxy 
stellar mass function from z — to z ~ 5. This model is therefore 
similarly well-suited to our study of fossil systems. If a halo in a 
merger tree has multiple progenitors, all but the most massive are 
considered to become subhaloes orbiting within the larger host halo 
and any galaxies they contain therefore become satellite galaxies in 
that halo. 

Due to the limited resolution of the Millennium Simulation 
(which may cause dynamical friction timescales to be poorly esti- 
mated), the time between becoming a satellite and merging with the 
central galaxy of the halo is computed from the analytic dynami- 
cal friction timescale. Specifically, each new satellite is randomly 
assigned orbital parameters from a distribution measured from N- 
body simulations and the appropriat e dynamical frictio n timescale 
computed following the approach of ICole et alj |2000), but multi- 
plied by a factor of 1. 5. This was found to produce the best fit to the 
luminosity function in Bower et al. ( 2006) b ut is also in good agree- 
ment with the results of Boylan-Kolchin, Ma & Ouataert ( 2008) 
who compared the predictions of analytic dynamical friction 
timescales with those from idealised N-body simulations. The 
satellite is allowed to orbit for the period of time calculated above, 
after which it is merged with the central galaxy of the halo. If the 
host halo doubles in mass before a satellite can merge, the satellite 
orbit is assumed to be reset by the merging which lead to that mass 
growth and so a new set of orbital parameters are assigned and a 
new merging timescale computed. 

2.3 The Millennium Gas Simulations 

The Millennium Gas Simulations are a series of hydrodynamical 
models constructed within the same volume, and values of initial 
perturbation amplitudes and phas es, as the parent dar k-matter-only 
Millennium Simulation (see, e.g. jHarflev et al.ll2008l) . Of the three 
principal models completed in this work, each contains additional 
baryonic physics: (i) the first does not follow the effects of radiative 
cooling and so overpredicts the luminosities of group-scale objects 
significantly, (ii) the second includes a simple preheating scheme 
that is tuned to match the observed X-ray properties of clusters at 
the present day and (iii) the third includes a simple feedback model 
that matches the observed properties of clusters today. We have 
used the second of these models in this work, as we only utilise 
the hydrodynamical properties of the groups at z = 0, where the 
observational and simulation results are well matched. 

The Millennium Gas Simulations consist of 5 x 10 s particles 
of each species, resulting in a dark matter mass of 1.422 x 10 10 /i _1 
Mq per particle and a gas mass of 3.12 x 10 9 /i _1 M Q per parti- 
cle. The Millennium Simulation has roughly 20 times better mass 
resolution than this and so some perturbation of the dark mat- 
ter halo locations is to be expected. In practice the position and 
mass of dark matter haloes above I0 13 h~ 1 M© are recovered to 
within 50 h' 1 kpc between the two volumes, allowing straightfor- 
ward halo-halo matching in the large majority of cases. 

The Millennium gas simulations used exactly the same cosmo- 
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logical parameters as those of the dark matter simulations. With the 
inclusion of a gaseous component, additional care needs to be taken 
in choosing th e gravitational softening l ength in order to avoid spu- 
rious heating dSteinmetz & W hite 1997). We use a comoving value 
of 2 5(l + z)/i~ kpc, ro ughly 4% of the mean inter-particle separa- 
tion teorgani et al.l2006l) until 2 = 3, above which a maximum co- 
moving value of 100/i -1 kpc is adopted. A different output strategy 
is followed in the Millennium Gas Simulations, where the results 
are output uniformly in time with an interval roughly correspond- 
ing to the dynamical time of objects of interest. This strategy results 
in 160 rather than 64 outputs and places particular emphasis on the 
late stages of the simulation. 



3 DATASETS USED IN THIS WORK 

We start with a catalogue of groups extracted by the friends-of- 
friends (FoF) algorithm employed in the Millennium Dark matter 
runs. Hereafter a "group" or "group halo" would refer to a group 
taken from this catalogue. In order to follow the evolution of groups 
from z~l to the present epoch, we have to combine various sets of 
information from this FoF group catalogue and the associated semi- 
analytic catalogues of galaxies, as well as the gas simulations. 

We select all groups of M(R 20 o) 5? 10 13 h' 1 M from the 
FoF group catalogue at z = 0.998. The mass cut-off is intended 
to ensure that the progenitors of the present day galaxy groups are 
indeed groups at z ~1 with at least four or five members (galaxies), 
above the magnitude cut of the catalogue. The evolution of each 
group was followed from z = 0.998 to z — (at 23 discrete 
values, equally spaced in log z) by matching the position of each 
halo to its descendants at later redshifts. 

The position of the central galaxy of each galaxy group, and 
the corresponding dark matter halo, were used to identify the mem- 
ber galaxies of each group. At each redshift and for each group 
halo, optical properties were extracted for its corresponding galax- 
ies from the semi-analytic galaxy catalogue. The model galaxies 
become incomplete below a magnitude limit of Mk — 5 log(/i) ~ 

— 19.7, due to the limited mass resolution of the Millennium simu- 
lation. We applied a /f-band absolute magnitude cut-off of Mk < 

— 19 on galaxies at all redshifts. During the matching process, for 
more than 99% of the groups at each redshift, corresponding galax- 
ies were found in the semi-analytic galaxy catalogue. The remain- 
ing groups were excluded from our final compiled list. 

In order to find the gas properties of all groups at z — 0, we 
cross-correlated our list of groups with the Millennium gas cat- 
alogue, and find the bolometric X-ray luminosity of our selected 
groups at 2=0. 

Out of 19066 dark matter group haloes with Af(7?20o) ^ 
10 13 h~ 1 M@ selected at z — 0.998, optical properties from the 
semi-analytic catalogue (as well as gas properties from the gas sim- 
ulations at 2=0) and the entire history of evolution at all redshifts 
up to 2=0, were found for 17866 (~94% of the initial sample at 
z ~ 1) of group haloes. Fig. Q] shows the bolometric X-ray lu- 
minosity from the Millennium gas simulation, plotted against the 
corresponding dark matter halo mass of each group, at redshift 2=0 
for all of the matched 17866 groups. 

The vertical dashed line in Fig. Q]corresponds to the conven- 
tional X-ray luminosit y threshold (Lx.b oi = 0.25 x 10 42 h~ 2 erg 
s~ ) for fossil groups dlones et al.|[2003h . as adopted in Sec. 14.2.11 
to define X-ray bright groups. There are 14628 groups above this 
threshold, out of 17866 groups. These X-ray bright groups will con- 
stitute the main data set for the rest of our analysis, except for £|4. 1 1 
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Figure 1. The relation between the mass of group haloes (within -R200) 
at z = from the Millennium DM simulation, and the bolometric X- 
ray luminosity of the corresponding haloes in the Millennium gas simu- 
lation. All groups have M(i? 2 oo) > !0 13 h^ 1 M Q at z ~ 1.0. The ver- 
tical dashed-line corresponds to the X-ray luminosity threshold Lx bol = 
0.25 X 10 42 h~ 2 erg s _1 generally adopted to define fossil groups (see 
Sec. |4~2~T1 . Of the 17866 groups matched in the two catalogues, 14628 
groups lie above this threshold. In this paper, we call these "X-ray bright 
groups". 



and ^5.31 where the whole range of halo mass will be explored 
to study the magnitude gap statistics and the local environment of 
groups, 



4 RESULTS 

4.1 The i?-band Magnitude Gap Statistic 

The dynamical friction /d yll will cause the more luminous galax- 
ies in a group to merge on a time scale which depends upon the 
velocity dispersion of the group, and since /d yn oc v~ 2 this is 
more freque nt in poorer groups than in clusters (e.g. iMiles et all 
l2004l 2006). As a result, on group scales, the likelihood of a few 
of the brightest galaxies merging to form the brightest galaxy, lead- 
ing to a large magnitude gap within a Hubble time, is higher. Thus, 
the distribution of the magnitude gap between the brightest galaxy, 
and the second and third brightest galaxies, in each group, is of- 
ten used as an indicator the dynamical age of group, particularly in 
fossil group s (Milosavljevic et al. 2006; van den Bosch "et"al .120071 : 
iDariush etal]|2007r . lvon Benda-Beckmann et al.ll2008h . 

We determine the magnitud e gaps from the Millennium semi - 
analytic models of iBower et all j2006t) and ICroton et all J2006h . 
and compare them with observational results from th e Sloan Dig- 
ital Sk y Survey (SDSS) C4 cluster catalogue data of Mill er et all 
j2005l) and the 2- degree Field Galaxy Redshi ft Survey (2dFGRS) 
group catalogue of Ivan den Bosch etal ]{2007l). 

The 2dFGRS group catalogu e is constru c ted ba sed on a halo- 
based group finder algorithm of lYang et alj d2005h and contains 
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~ 6300 groups within the mass range log(Af//i _1 Mq) ^ 13.0 
where the mass of each group has been determined from the total 
luminosity of all group members brig hter than Mtj — 5 log h = 
-18. The C4 catalogue iMiller et alj ( f2005l) consists of ~ 730 
clusters identified in the spectroscopic sample of the Second Data 
Release (DR2) of the SDSS inside the mass range 13.69 ^ 
log(A///i _1 Mq) ^ 15.0, estimated from the total r-band optical 
luminosity of cluster galaxies. 

The results of the comparison between our distribution of 
the es timated R-hmA ma gnitude gap s from semi-analyti c mod- 
els of iBoweret aD (§006) as well as ICroton et alj i feOOd) . based 
on the Millennium simulation (red triangles and black circles re- 
spectively), and the observed results from the C4 cluster catalogue 
and the 2dFGRS group catalogues (blue histogram), are shown 
in Fig. [2] The ma gnitude gap statistics Ami2 and Amu from 
lBoweretalJj2006l) are in excellent agreement with those obtained 
from 2dFGRS group catalogue and SDSS C4 catalo gue of clus- 
ters. H owever, the semi-analytic galaxy catalogue of ICroton et al.l 
(2006) predicts a larger fraction of groups with Am^ ^ 2.0 
for both the SDSS and 2dFGRS samples. This is in particular of 
great importance to the determination of the space density of fossil 
galaxy groups, and the comparison of fossil samples drawn from 
simulated and observed catalogues, whic h use the magnitude gap 
as a key discriminant, (e.g., see Table 1 of lDariush etaU 20071) . 

The shift of the distrib ution of Amu and Amy to larger 
values than observed in the Croton et al. d2006h model may re- 
flect the fact that in the lBower et alj ( 12006]) m odel, the treatment o f 
dynamical friction differs from that used by ICroton et al. I d2006h . 
In both models, N-body dynamics are used to follow the orbital 
decay of satellite galaxies whose subhalo can be resolved. How- 
ever, when the subhalo can no longer be re liably followed, the 
dynamical friction calculations differ. In the iBower et al. I d2006t) 
model, the dynamical f riction timescale is initially calculated fol- 
lowing |CoieetaE| d200d) . However, if the host halo of the satellite 
is deemed to undergo a formation event (corresponding to a mass 
doubling since the previous formation event), before the satellite 
merges, then a new orbit for the satellite is selected at random, 
and the dynamical friction timescale for the satellite is recalcu- 
lated. This calculation takes into account the scattering of galaxies 
to larger energy orbits during the merger of their pare nt halo. 

Another possible cause for the success of the IBower et al.l 

d2006l) model in matching the magnitude gap statistics is that it 
predicts a large scatter in the relation betwee n galaxy stellar mass 
and halo mass — significantly more than in the lCroton et ai] d2006l) 
model. As a result, sometimes rather large satellite haloes arrive 
carrying relatively small galaxies resulting in a large difference be- 
tween the magnitude of the dominant object and the next most lu- 
minous. This difference occurs because the AGN feedback is not 
gu aranteed to switch o ff the cooling at a particular halo mass in 
the IBower et alj (2006) model, as it depends on the merging and 
cooling history of each halo. 

For the purposes of this work, the fact that the IBower et alj 

(2006) method for computing merging timescales results in good 
agreement with the observed magnitude gap distributions, makes it 
well suited for the remainder of our study. 

4.2 Evolution of galaxy groups 

In cosmological simulations, the age of galaxy groups can be ex- 
pressed in terms of the rate of the mass assembly of the groups. This 
means that for a given group halo mass, groups that formed early, 
assemble most of their masses at an earlier epoch in comparison to 



younger groups. Thus the assembly time of a dark matter halo, de- 
fined as the look-back time at which its main progenitor reaches a 
mass that is half of the halo's present day mass, is larger in "older" 
systems than in "younger" ones. Of course, in cosmological simula- 
tions such as the Millennium runs, where the structures in Universe 
from hierarchically, massive systems which form later turn out to 
have shorter assembly time than low mass groups. Therefore one 
should take into account the mass of systems when comparing the 
mass assembly of various types of groups and clusters. 

4.2.1 Fossil groups of galaxies 

How is the history of mass assembly of a group or cluster related to 
its present observable parameters? It is expected that groups which 
have formed earlier tend to be more dynamically relaxed, thus re- 
sulting in a hotter intergala ctic medium (IGM) and being mor e 
likely to be X-ray luminous dMiles et alj2004lForbes et alj2 006). 
It has been shown that in X-ray luminous systems, the brightest 
galaxies are bigger and more optically luminous and those be- 
longin g to systems that have little or no diffuse X-ray emission 
(e.g.. lKhosroshahi et alj|2004l) . On the other hand, groups with the 
same mass that have formed late, and are still in a state of col- 
lapse, would not show X-ray emission associated with their IGM 
(Rasmussen et al. 2006; lBaietai] [2010), and are less likely to be 
dominated by a massive elliptical in their cores. 

Hitherto the so-called fossil galaxy groups, which are sup- 
posed to be canonical examples of groups that have formed 
early, have been identified by requiring that their X-ray lu- 
minosities exceed Lx.boi ^ 0-25 x 10 42 /t~ 2 erg s -1 (e.g. 
iKhosro shahi. Ponman, & Jones 2 0071 ; Ijones et al] |2003l) . In addi- 
tion, a fossil group needs to have, within half a virial radius of 
the group's centre, the second brightest galaxy to be at least 2- 
mag fainter than the brightest galaxy, i.e. Amu ^ 2.0 Q. So far 
these two observational criteria have been jointly used to explore 
fossil groups and clusters of galaxies. Therefore, Fig. [3] which dis- 
plays all the X-ray bright groups fulfilling the X-ray criterion in 
Fig. Q] an d the optical criterion Amu ^ 2.0 (dotted horizontal 
line) should separate groups which have been formed earlier in 
comparison to their counterparts with Amu < 2.0. 

Note that, in numerical simulations, fossils are identified as 
groups with Amu *s 2.0 within R200 or 0.5 Juno- Our results 
from this study as well as those represented in lDariush et alj d2007l) 
show that the fraction of fossils (and therefore their space densi- 
ties) depend on the search radius within which Amu is estimated, 
whereas the history or mass assembly does not change that much. 

4.2.2 The mass assembly of X-ray fossil groups 

Let us introduce the parameter a z which for an individual group 
is the ratio of its mass at redshift z to its final mass at z = 0, i.e. 
a z = M z /M z -o- Thus at a given redshift z, groups with larger 
a z have assembled a larger fraction of their final mass by z than 
groups with smaller values of a z . 

In Fig.[3] we plot the magnitude gap Amu (within 0.5iuoo), 
estimated for all 14628 X-ray bright groups (i.e. groups with 
£x,bol 0.25 x 10 42 /i _2 erg s _1 ) at z — as a function of their 
mass fraction ai.o at z = 1. Groups are colour-coded according to 
their dark matter halo mass. The horizontal dashed line separates 

1 This condition can be replaced by log(L2/Li) ^ —0.8 where L\ and 
L2 are the luminosities of the first two brightest galaxies. 
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Figure 2. The ij-band magnitude gap distribution for haloes from the Millenni um semi-analytic models oflBower et all 120061) (red triangles) and Crotonetal 
1 2006) (black circles) superposed on the data from 2dFGRS group catalogue o flvan den Bosch etail 120071) as well as SDSS C4 cluster catalogue o fjMiller et all 
12005) (blue histograms), (a) The magnitude gap A77112 between the the first and second most luminous galaxies, compared with galaxies from the SDSS 
C4 catalogue of clusters computed within projected radius of 5Q0h~ x kpc. (b) The same as in (a) but for the magnitude gap Am^ between the first and the 
third most luminous galaxies, (c) The magnitude gap Amu estimated within -R200. compared with galaxies from the 2dFGRS group catalogue. The ~ 6300 
2dFGRS groups are within the mass range \og(M (R200) / h~ 1 Mq) ^ 13.0, and those from SDSS C4 catalogue consist of ~ 730 clusters within mass 
range 13.69 ^ log(M(i? 2 oo)M _1 M ) s£ 15.0. 



groups into fossils (Amu ^ 2.0) and non-fossils (Amu < 2.0). 
All data points on the right side of the vertical dashed-line have as- 
sembled more than 50% of their mass by z ~1.0 and hence have a 
minimum assembly time of about ~7.7 Gyr. The contour lines rep- 
resent the number of data points (groups) in each of 25 x 25 cells of 
an overlaid grid which is equally spaced along both the horizontal 
and vertical axes. 

Three results emerge from this plot: (i) As is expected, on av- 
erage the rate of mass growth in massive systems is higher than 
in low mass groups as the majority of massive groups and clusters 
have assembled less than 50% of their final mass at z ~1.0. (ii) 
Less massive groups (and therefore early-formed ones) tend to de- 
velop larger magnitude gaps in comparison to massive groups and 
clusters. Consequently the fraction of massive fossils, identified in 
this way, is less than low mass fossil groups, (iii) For any given 
Qi.o <; 0.5, the majority of groups have magnitude gaps Amu < 
2.0, as is evident from the density of contours. In other words, the 
number of early-formed groups with values of Amu < 2.0 ex- 
ceeds the number of fossil groups with Amu > 2.0. 

Unlike the first two results, the third conclusion is not in agree- 
ment with our current view that early-formed groups necessarily 
develop larger magnitude gaps. Clearly, the majority of groups with 
similar values of ai.n > 0.5 have smaller magnitude gaps. Without 
doubt, the parameter Am 12 is influenced by the infall and merg- 
ing of galaxies and sub-groups within galaxy groups. This could 
result in the increase (in case of merging) or decrease (in case 
of the infalling of new galaxi e s) in Amu- Indeed, in the work 
of Ivon Benda-Beckmann et alj j2008), one finds that the "fossil" 
phase of any fossil group is transient, since the magnitude gap cri- 
terion will sooner or later be violated by a galaxy comparable to the 
brightest galaxy falling into the core of the group. 



To quantify the above results, we study the evolution with red- 
shift of various physical parameters for two different sample of 
groups, drawn from the distribution of galaxy groups in Fig. [3] In 
the first sample, haloes are divided into old(oi.o 5^ 0.5, b + c in 
Fig. [3) and young (ai.o ^ 0.5, a + d in Fig. [3) groups respec- 
tively. In the second population, haloes are classified as X-ray fos- 
sil (Amu ^ 2.0, a + b in Fig. [3} and control (Amu ^ 0.5, c + d 
in Fig. [3} groups based on the magnitude gap between the first and 
the second brightest galaxies within half a virial radius of the centre 
of the group. 

For each sample, the evolution of various parameters are 
shown in two panels of Fig. [4] From top to bottom these param- 
eters are: 

• a z , i.e. the ratio of the group halo mass at redshift 2 to its mass 
at z = 0, 

• Ami2 within 0.5i?200, 

• Ratio of the number of galaxies within 0.5i? 2 oo at redshift 2 
of a given galaxy group to the number of galaxies within 0.5J?200 
at redshift z = of the same group, i.e. G z , 

• Group velocity dispersion ay in km s~ . 

In each panel of Fig. [4] the left, middle, and right columns 
correspond to different ranges in group mass, as indicated. The left 
panel in Fig. [4] illustrates X-ray fossil (red triangles) and control 
(blue circles) groups respectively while the right panel shows old 
(red stars) and young (blue squares) groups. The horizontal green 
dashed-lines intersect y-axes at G z = 1. Errors on data points are 
the standard error on the mean, i.e. a/ %/~N, where a is the standard 
deviation of the original distribution and N is the sample size. 

A comparison between Figs.|4jil, a2, a3 and Figs.|4ji4, o5, a6 
shows that older groups, which have been picked up according to 
their lower rate of mass growth (i.e. larger ai.o), represent a per- 
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Figure 3. The magnitude gap Ami2 within 0.5 -R200. estimated for all 14628 X-ray bright groups in Fig.[T](i.e. groups with Lx.bol ^ 0.25 X 10 42 /i — 2 erg 
s _1 ) at z = versus the ratio of the group halo mass at redshift z = 1 to its mass at z = (ai.o). The horizontal dashed-line separates groups into fossils 
(Amj2 ^ 2.0) and non-fossils (Amu < 2.0). The vertical dashed-line corresponds to ai.o = 0.5. Groups with ai.o 0.5 have formed more than half 
of their mass by z ~ 1.0 and hence have a minimum assembly time of about ~ 7.7 Gyr. Data points are colour-coded according to FoF group halo mass 
Mr200 a t present epoch. The density of data points is represented by black contour lines which is the number of groups in each of 25 X 25 cells of an overlaid 
grid, equally spaced horizontally and vertically. The upper panel represents the histogram of X-ray bright fossil groups, i.e. all groups with Amu ^2.0 and 
ix.bol > 0.25 X 10 42 /i- 2 ergs- 1 . 



feet class of fossils, though they develop a magnitude gap Ami2 
which is not as large as those seen in X-ray fossils (see also 
Figs.|!>l,...,&6). 

On the other hand, unlike old groups, X-ray fossils de- 
velop large magnitude gaps, which do not necessarily correspond 
to their rapid mass growth, especially in massive groups with 
log(M(R 2 oo)/h~ 1 M Q ) ^ 14.0. This reflects the fact that the ma- 



jority of real passive groups have a small magnitude gap between 
their two brightest galaxies. Thus the expression Amu ^ 2 only 
partially separates genuine old/passive groups from young/forming 
groups, as there are a larger fraction of genuine old groups but with 
small Ami2. 

From Figs.|4j;4, c5, c6, it is obvious that older groups are es- 
sentially more relaxed systems, which have not recently experi- 
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Figure 4. The evolution with redshift of various physical parameters of X-ray bright groups, in various ranges of group mass. Left panel: Haloes are classified 
as X-ray fossil (Amu 5S 2.0, red triangles) and control (Amu ^ 0.5, blue circles) groups based on the magnitude gap between the first and the second 
brightest galaxies within 0.5i?200- Right panel: Groups are divided into old(ai.o ^ 0.5, red stars) and young (eti.o ^ 0.5, blue squares) population 
respectively. Each row represents the evolution of a parameter characteristic of galaxy groups that can be measured from the simulations (but not necessarily 
from observations). From top to bottom the parameters on the y-axis represent: (i) a z , the ratio of the group halo mass at redshift z to its mass at z = 
(al, ...,a6), (ii) Amu, i.e. the magnitude gap between the first two brightest group galaxies found within 0.5/?200 ...,b6), (iii) G z , the ratio of the 
number of galaxies within 0.5 -R200 at redshift 2 of a given galaxy group to the number of galaxies within 0.5 -R200 at redshift 2 = of the same group 
(cl, ...,c6), and (iv) the group velocity dispersion ay in km s _1 (dl, ...,d6). In the third panels from top, the horizontal green dashed-Iines intersect the 
y-axis at G z = 1. 



enced a major merger, as the rate of infall of galaxies is equal 
or even less than the rate at which galaxies merge with the cen- 
tral group galaxy. Therefore in old groups the parameter G z is 
more or less constant with time, compared to that of the younger 
groups within the same group mass bin. The situation is a bit 
different in X-ray fossils with log(A/(i?2oo)/ft _1 M ) sC 14.0 
(Figs. Hpl, d2), since, in these cases, the rate of galaxy merging 
is noticeably larger than that of the infall of galaxies. As a result, 
very large magnitude gaps can develop in X-ray fossil groups. It 
is also evident both massive X-ray fossils and control groups with 
log(M(i?2oo)/ft~ 1 M ) ^ 14.0 (Fig.gJS) are in a state of rapid 
mass growth. As a consequence, massive X-ray groups are not dy- 
namically relaxed systems as they are influenced by infall of galax- 
ies and substructures. 



Finally, it is worth considering how the velocity dispersion, 
plotted in Figs.|4jil, d6 changes with time in the different kinds 
of groups. As Figs.|4]i4, d5, d6 show, as long as the rate of infall of 
subgroups is close to 1.0 (green dashed-line), the velocity disper- 
sion does not significantly change with time, which in turn is a sign 
that these groups are certainly relaxed systems. 



4.2.3 BCG magnitudes 

Since the central galaxy in a fossil group is a product of numer- 
ous mergers, many of them with luminosities close to L* galaxies, 
X-ray fossils are expected to be dominated by optically luminous 
brightest galaxies (BCGs) more often than their non-fossil counter- 
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Figure 5. The absolute R-band magnitude of BCGs for all X-ray bright groups versus Amu within O.5R200 in four different mass bins. A grid of 45 X 55 
cells has been superposed in each panel, and the cell is colour-coded according to the median value of a 1.0 within the cell. The contours trace the distribution 
of median ai.o values in each panel. All panels have the same scale. 



parts. Here we explore the correlation between the luminosity of 
the central galaxies of groups with large magnitude gaps, and their 
mass assembly history. 

The four panels in Fig. [5] demonstrate the relation between 
the absolute .R-band magnitude of the BCGs for all X-ray bright 
groups and the magnitude gap Amu within 0.5 R200 in four dif- 
ferent mass bins. In each panel of Fig.[5] a grid of 45 x 55 cells is 
overlaid, where each cell is colour-coded according to the median 
of Q1.0 in that cell. Accordingly, contours in each panel trace the 
distribution of median ai.o values. 

Fig. [5] shows clearly that both ai.o and BCG .R-band magni- 
tudes increase with decrease in group halo masses. But it does not 
show a tight correlation between the R-band luminosity of BCGs 
and group magnitude gaps Amu, though the correlation is more 
pronounced in clusters with M{R 2 qo) > 10 14 h' 1 M (Fig.[5ji). 
Therefore, imposing a magnitude cut for the BCGs would result in 
the loss from the sample of a large number of genuine groups which 
are not X-ray fossil systems according to the optical condition in- 
volving Amu. 



4.3 The Fossil phase in the life of groups 

The existence of large magnitude gaps in X-ray fossils in 
Figs.[4j>l, 62, 63 is expected as these groups were initially selected 
according to their Amu at z=0. It would be interesting if they 
could be shown to have maintained such large magnitude gaps for 



a longer time, in comparison to control groups, which would be 
the case if X-ray fossils were relaxed groups without recent ma- 
jor mergers. Also if fossil groups in general are the end results of 
galaxy merging, then we do expect the majority of fossils selected 
at higher redshifts to still be detected as fossils at the present epoch. 
In other words, the fossil phase in the life of a galaxy group should 
be long-lasting. 

To test this, we select three sets of fossil groups with Amu 
2.0 at three different redshifts. By tracing the magnitude gap of 
each set from z=1.0 to z=0, we examine the fossil phase of each 
set in time. In Fig. [6] fossils (black shaded histogram) and control 
(black thick line histogram) groups are selected at z=Q (left col- 
umn), 2=0.4 (middle column), and z=1.0 (right column). Fractions 
of fossil and control groups in each column of Fig. [6] have been 
separately estimated by normalising the number of fossil and con- 
trol groups at other redshifts to their total numbers at the redshift at 
which they were initially selected. 

This plot shows that, contrary to expectation, no matter at what 
redshift the fossils are selected, after ~4 Gyr, more than ~90% of 
them change their status and become non-fossils according to the 
magnitude gap criterion. Over the span of 7.7 Gyr, which is the time 
interval between z =0-1, very few groups retain a two-magnitude 
gap between the two brightest galaxies. This means that the fossil- 
phase is a temporary phase in th e life of fossil groups (also see 
Ivon B enda-Beckman n et al. 2008), and there is no guarantee that 
an observed fossil group, at a relatively high redshift, remains a 
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Figure 6. The fate of fossil groups identified at different redshifts. Fos- 
sil groups (dark shaded histogram) and control groups (grey shaded his- 
togram) are initially selected at 2=0 (left column), 2=0.4 (middle column), 
and 2=1.0 (right column). For these objects, we explore what fraction re- 
main "fossil" or "control" groups at two other epochs in the redshift range 
< 2 < 1.0. It is clear that the fossil phase does not last in > 90% of 
groups after 4 Gyr, no matter at which epoch they are identified. 



fossil until the present time, if fossils are selected according to their 
magnitude gap Ami2 > 2.0. 



5 REVISING THE OPTICAL CRITERION FOR FINDING 
FOSSIL GROUPS 

Using the Millennium simulation DM runs as well as the gas and 
semi-analytic galaxy catalogues based on them, it seems from the 
above that the conventional optical condition Amy ^ 2, used to 
classify groups as fossils, does not ensure that that these systems 
represent a class of old galaxy groups, in which the central galaxy 
has grown through the merging of other comparable group galaxy 
members. Having said that, it is true that the magnitude gap in a 
galaxy group is related to the mass assembly history of the group, 
for we saw that in groups, such a gap develops gradually with time. 
However, the difference between the luminosities of the two bright- 
est galaxies in groups is not always reliable for the identification of 
fossil systems, as this quantity is vulnerable to the assimilation of 
a comparable galaxy into the core of the group, as a result of infall 
or merger with another group. 

We therefore attempt to identify a more robust criterion, in 
terms of the difference of optical magnitudes among the brightest 
galaxies in a group, which might be better suited to identifying sys- 
tems where most of the mass has been assembled at an early epoch. 
As introduced in Sec. 14.2.21 we quantify the age of a group in terms 
of the mass assembly parameter ai.o, which for an individual group 
is the ratio of its mass at redshift z — 1.0 to its present mass at z — 0, 
i.e. ai.o = M z -i.o/M z -o- We begin by considering the effect of 
the radius within which the magnitude gap is calculated. 



5.1 A general criterion for the magnitude gap 

Assume a general optical condition in defining early-formed groups 
according to the magnitude gap between the brightest group galaxy 
and other group members in the following form: 



Amu > j, 



(1) 



where Amu is the difference in i?-band magnitude between the 
first brightest group galaxy and the i th brightest group galaxy 
within 0.5R.200 (or R200) of the group centre. The current definition 
of fossils involves i — 2 and j = 2. Obviously any group satisfying 
Eq. Q]must contain at least i galaxies. We do not consider i > 10 
since then we have to exclude most groups in our sample, and it 
would turn out not be very useful for observers as well. 

As we consider the magnitude gap between the brightest to 
the i th (=2, 3,..., 9, 10) brightest galaxy, the value of the magnitude 
gap varies from j >0 to j <5. Our aim is to find a pair of (i, j) in 
Eq.Q] which yields the best selection of genuinely old groups with 
a history of early mass assembly. 

In Fig. [7] we show how the parameter ai.o = M z= i.o/M z= o, 
which represents the mass assembly of groups since redshift z=1.0, 
depends upon the selection of i and j in Eq.Q] For each value of i 
in Fig. [7] groups are first sorted according to their magnitude gaps 
Amu estimated within O.5R200 or R200), where the latter is the 
overdensity radius of the group. For each i, the average value for 
ai.o is calculated, for each value of j, for all groups satisfying 
Eq.Q] The plot is colour-coded according to the values of ai.o. 
The black contours give an idea of the number of groups involved: 
the fraction of the total number of groups identified by parameters 
(i, j) is constant along each of these lines. 

From Fig. |7Ji, for instance, we find that systems with i — 4 
and j — 3, i.e. systems for which Amu *5 3 (within O.5R200), 
yield ~2.4% of groups with ai.o ~0.54. The same is ~1.2% with 
ai.o ~0.56, if one estimates Amu ^ 3 within R200 according to 
Fig. |7p. 

In fact, by changing our search radius from O.5R200 to R200, 
we find ~50% fewer groups satisfying Eq.Q] for the same value of 
ai.o. Therefore, hereafter, we just use Panel a of Fig. [7j which es- 
timates the magnitude gap Amu within half of the overdensity ra- 
dius, i.e. 0.5R200- From this plot, the fraction of groups picked out 
by applying the conventional fossil criterion Amu ^ 2 is ~4.0% 
with ai.o ~0.52. 

If we were to find an improved criterion for finding fossils, a 
better set of parameters (i, j) in Eq.Q]should 

(i) identify groups with larger valuer of ai.o, and/or 

(ii) find a larger fraction of groups with the same or larger value 

of ai.o, 

than found in conventional fossils, i.e. groups with (i, j)=(2, 2). 

For example, by choosing i — 4 and j = 3, the fraction of 
groups found with Amy ^ 3 turns out to be ~40% less than 
when i=j=2, but it would identify slightly older groups, with an 
average ai.o ~0.55, whereas the average ai.o ~0.52 in fossils 
with i=j=2. In other words, (i, j)=(4, 3) identifies marginally older 
groups, at the expense of losing a large number of early-formed 
groups, compared to the case of conventional fossils (i, j)=(2, 2). 

Exploring Fig. |7Ji, we adopt (i, j)=(4, 2.5) as an example of 
how the fossil search criterion can be improved. If we define all 
groups with Amy 2.5 within 0.5 -R200 as fossils, then the we 
would find groups with on average the same mass assembly history, 
i.e, the same average value of ai.o, but we would identify ~50% 
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Figure 7. The dependence on (i, j), defined in Eq.fJ] of the mass assembly parameter cti.o, which is defined as the ratio of the mass of a group at redshift 
z = 1.0 to its mass at z = 0, i.e. ai.o = M z= i.q/M z =0' For each value of i, groups are sorted according to the value of their magnitude gaps Aroii 
calculated within a certain radius (different for the two panels). Then, for different values of j, the average for ai.o is calculated for all groups satisfying Eq.fJ] 
The plot is colour-coded according to ai.o. The black contours are drawn such that the fraction of the total number of groups identified is constant along each 
line. The magnitude gap Amu is calculated (Panel a.) within half the overdensity radius, i.e. 0.5R-2OO, and (Panel b.) within the overdensity radius, i.e. R200- 



more such groups, compared to groups identified the conventional 
parameters (i, j)=(2, 2). 

We will examine such groups further in the next section. 
Meanwhile, the Figs. [7] will allow the user to find their favourite 
combination of (i, j) for both O.5R200 and R.2oo- 



5.2 The optical criterion Am^ ^2.5 within 0.5R 200 

Having explored alternative criteria for identifying groups with a 
history of early formation, we now compare the history of mass as- 
sembly of groups selected according to A77114 ^ 2.5 (these groups 
are hereafter collectively referred to as F14) with those groups 
selected according to Am^ 2.0 (hereafter F12), both within 
O.5R200 of the group centre. The latter category are the conven- 
tional fossil groups. 

The blue histogram in Fig.[8]represents the distribution of the 
mass assembly parameter ai.o for all X-ray bright groups (as de- 
fined in Fig.QJ in our sample. It also shows the groups in the cate- 
gories F14 {red thick histogram) and F12 (grey shaded histogram). 
The green shaded histogram corresponds to groups that satisfy both 
criteria, i.e. F12 n F14. The green dash dotted line intersects the 
x-axis at ai.o = 0.5, representing groups for which half of their 
mass had been assembled at redshift z = 1. Gaussian fits to each 
histogram are overlaid. 

Panels a, b, and c in Fig. [8] correspond to different ranges 
of the logarithm of the group halo mass 13.0 ^ logM(i?20o), 
13.0 < logM (i?20o) < 13.5, and logM (#200) ^ 13.5 respec- 
tively, where M (-R200) is in units of hT 1 Mq. These figures, as 
well as the values of the peaks of Gaussian fits to the distribution 
of ai.o in each case (given in TableQJ lead to the following obser- 
vations: 

(i) The groups belonging to F12 and F14 are older than the over- 



all population of X-ray bright groups (for all values of halo mass), 
though such a difference is less pronounced in low-mass systems. 
Since haloes are thought to be hierarchically assembled, one ex- 
pects to find a higher incidences of early-formed and low-mass 
groups in comparison to massive systems. 

(ii) Within the errors, groups belonging to F14 are almost as old 
as F12, i.e. the estimated ai.o in F12 systems (as given in Table. 
Q} is more or less the same as found in F14 systems. However, the 
fraction of groups in category F14 is at least 50% more than F12. 
This tells us that in general, the criterion Ami4 ^ 2.5 has a higher 
efficiency of identifying early-formed systems than A77112 ^ 2.0. 

(iii) Interestingly, ~ 75% of F12 haloes in Fig. [8^ also fulfil the 
Ami4 ^ 2.5 condition. Conversely, ~ 35% of F14 haloes satisfy 
the Ami2 ^ 2.0 criterion. This means that a large proportion of the 
population of early-formed groups in the category F14 is different 
from those in F12. Groups which satisfy both criteria, i.e. F12 D 
F14, are not necessarily older in comparison to those belonging to 
either F12 or F14 (see Table. [T). 

(iv) Fig.[8p shows that in fact neither criterion Amu ^ 2.0 nor 
Ami4 ^ 2.5 is efficient in finding early-formed groups in low- 
mass regime, even among X-ray bright groups. 

In the following section, we compare the environment and 
abundance of the groups belonging to the F12 and F14 categories. 



5.3 The local environment of fossil groups 

If galaxy mergers are responsible for the absence of bright galax- 
ies in groups such as X-ray bright fossils, then most of the matter 
infall into these systems would have happened at a relatively ear- 
lier epoch. Consequently, at the present time, old groups should 
be more isolated than groups which have recently formed (e.g. 
lLaBarbera et alj|2009h . Here we examine the local environment 
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Table 1. Peak values (from Gaussian fits) of the histograms of the mass assembly parameter ci\ o for 
various classes of groups, (see Fig. [8). Group halo mass M(i?20o) is in units of h -1 Mq. F14 consists of 
all groups with Amy ^ 2.5, and F12, those with Amu ^ 2.0, both within O.5R200 of the group centre. 



Group type 


13.0 log M(R 2 oo) 
Panel a 


13.0 s: logM(R 2 oo) «S 13.5 
Panel b 


log M(R 200 ) > 13.5 
Panel c 


All X-ray bright groups 


0.41 ±0.01 


0.52 ±0.01 


0.38 ± 0.01 


F12 


0.53 ±0.01 


0.56 ±0.01 


0.49 ± 0.01 


F14 


0.52 ±0.01 


0.55 ±0.01 


0.48 ± 0.01 


F12 n F14 


0.54 ±0.01 


0.55 ±0.01 


0.50 ±0.02 



log(M/M,h"')>13.0 



13.0<log(M/M r ,h" 1 )<13.5 
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Figure 8. Histograms of the mass assembly parameter cti .0 for X-ray bright groups (blue histogram), groups that satisfy the criterion F14 (red thick histogram), 
those that satisfy F12 (gray shaded histogram), and groups that satisfy both criteria, i.e. F12 n F14 (green shaded histogram). Overlaid are Gaussian fits to each 
histogram (see Table[TJ. Panels a, b, and c correspond to the logarithm of the group halo mass in the range 13.0 ^ log M(i?20o). 13.0 ^ log M(i?20o) ^ 
13.5, and log Af(i?20o) ^ 13.5 respectively, the unit of M(i?20o) being h~ 1 Mq. The green dash dotted line intersects the rr-axis at 01. = 0.5 where 
haloes have assembled 50% of their mass at redshift z=1.0. 



of groups, using the density parameter A4, defined as the number 
of haloes within a distance of 4 h^ 1 Mpc from the centre of each 
group. The the local densities are calculated at 2 = according to 



A 4 = ^-l, 

Pbg 



(2) 



A 4 



6.31 ±0.17 
6.25 ±0.18 
5.10 ±0.35 
5.19 ±0.26 



for 
for 
for 
for 



Control groups (Amu ^ 0.5) 

X-ray bright groups 

F12 y ' 

F14 



where p 4 is the number density of haloes within a spherical vol- 
ume of of 4h~ Mpc in radius, and pf, g is the background den- 
sity of haloes within the whole volume of the Millennium sim- 
ulation. Since the mass assembly of groups is mostly influenced 
by the infall of subgroups, which individually have masses typ- 
ically below ~10% (and often substantially smaller) of the par- 
ent halo mass, it is important to take into account all haloes with 
M (-R200) > 10 11 h" 1 M Q from the FoF group catalogue in order 
to estimate A 4 . 

From Gaussian fits to the histograms of the local density A 4 
of F12 and Fi 4 , control groups, as well as those of all X-ray bright 
groups, we find 



where A 4 is estimated using Eq. [2] It seems that both F12 and F14 
groups are more likely to lie in lower density regions than control 
groups and X-ray bright groups. This is in agreement with our ex- 
pectation as early-formed groups are assumed to be in low-density 
local environments. However, the local density around F12 and F14 
groups is more or less the same. This is consistent with the fact that 
F12 and Fi 4 have similar values of ai.o (see Table[TJ. 

A similar trend is seen for density measures A5 and Ag, but as 
the sampling volume increases (> Ae), the above trend disappears, 
showing that the trend is related to the immediate environment of 
groups. 
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5.4 The abundance of fossil groups 



Various studies have shown that the fraction of early-formed groups 
increases as the group halo mass decreases (e.g. lMilosavlievic et al.l 
12004 iDariush et af] |2007). This phenomenon reflects the fact that 
structures form hierarchically, where small virialised groups form 
early, whereas most massive clusters form late. As the merging of 
galaxies in clusters is less efficient than in groups, due to the high 
velocity dispersion of cluster galaxies, clusters are less likely to de- 
velop lar ge magnitude gaps. At t he same time, in low-mass groups 
(see.e.g jMiles et all2 004. 2006) dynamical friction is more effec- 
tive in ensuring galaxies fall to the core of the system, due to the 
smaller relative velocities involved. As a result, the existence of 
large magnitude gaps should be more frequently found in groups 
than in clusters. Thus, to find an old population of groups accord- 
ing to some criterion, and to study the way the criterion depends 
on group halo mass, would be a good test for the validity of the 
condition. 

The top panel of Fig. [9] displays the abundance of F12 (gray 
shaded histogram) and Fibred thick line) groups, defined as the 
fraction of haloes in each category, as a function of halo mass. The 
range of halo mass explored is logM(i?20o) > 13.4 in units of 
h^ 1 Mq. Below this mass limit, the number of groups abruptly de- 
crease, since all groups here have been chosen to be X-ray bright 
groups (see Fig.[TJl. 

The plot shows that in comparison to the F12 groups, the F14 
groups are populated by less massive haloes. This can be better 
seen in the bottom panel of Fig.[9] where the relative fraction of F14 
groups over F12 is shown. It can be inferred that, on average, the 
fraction of F14 groups with halo mass M(i?20o) =Ss 10 14 M© 
is at least 50% more than the fraction of F12 . However, in the 
mass range M(i?20o) ^ 10 14 M©, the fraction of F14 groups 
decreases, though since the overall numbers in the extremely high 
mass range are low, the statistics are poorer. 



5.5 The survival of the magnitude gap: F12 vs. F14 



In Sec. 14. 31 it was found that, in general, for the conventional fossil 
groups (F12), the /o.s.«7 phase is transient, > 90% of such groups 
ceasing to remain fossils after 4 Gyr. Here we examine whether the 
fossil phase in F\& groups fares better. 

The histograms in Fig.[lO]represent the fractions of F12 (black 
line) and F14 (thick red line), as a function of look-back time in 
Gyr. The plot shows that in comparison to the F12 groups, fhe/as- 
sil phase lasts longer by almost 1 Gyr for the same fraction of F14 
groups. For example, the fraction of F12 groups that maintains its 
magnitude gap after ~2.2 Gyr, falls to 28%, while the correspond- 
ing period is ~3.2 Gyr in F14. Thus, not only does the Amu J? 2.5 
condition identify at least 1.5 times as many fossil groups as the 
F12 condition, it also identifies groups in which the fossil phase 
lasts significantly longer. This can be explained from our analysis 
of the halo mass distribution within F12 and F14, already discussed 
in Sec. E4] 



5.6 Comparison with observed groups 

When making detailed comparisons between simulations and cata- 
logues of galaxies and groups compiled from observations, one has 
to be aware that simulated dark matter haloes have limited resolu- 
tion, as mentioned in § 12.2. H and 12.2.21 This means that all galaxies 
in the semi-analytic models, assigned to a particular halo, might not 
belong to dark matter subhalos. We find that even after applying a 
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Figure 9. Top panel: The abundance of Fi2(gray shaded histogram) and 
Fibred thick line) groups, i.e., the fraction of groups in each category as a 
function of halo mass. Bottom panel: The relative fraction of F14 over F12 
groups as a function of halo mass. 



magnitude cut, there would be a significant number of modelled 
galaxies, whose orbits are analytically calculated, that would end 
up being not a member of a sub-halo, and thus would not be classi- 
fied as a group member. 

While dealing with magnitude gaps Amu and Amu, it is 
worth examining to what extent these quantities are vulnerable to 
resolution effects like the above. A direct way would be to com- 
pare the magnitude gap distribution of groups, selected based on 
their Amu, between simulations and observations. This is not a 
straightforward task, as groups identified from observational sky 
surveys are biased due to incompleteness in measured magnitude 
and redshift. Furthermore, a variety of group finding algorithms 
are been adopted to identify groups in simulations and observation, 
which adds uncertainties to any such co mparison. 

Here, we use the group catalogue o f lYangetaljj2007l) , which 
uses a halo-based group finder on the Sloan Digital Sky Survey 
(SDSS DR4). They define groups as systems whose dark matter 
haloes, have an overdensity of 180, determined from dynamics. 
This makes this catalogue suitable for comparison with the Mil- 
lennium simulation, where dark matter halos have an overdensity 
of 200. From Sample II of the catalogue, groups with following 
properties are selected: 

• they have at least four members, 

• they are within the redshift range 0.01 ^ z ^ 0.1, and 

• their estimated halo mass is log M (.Ri 80 ) > 13.25, h- 1 M Q , 
since our Millennium X-ray groups have a similar mass threshold 
(see Fig.[T}. 

After applying the above criteria, 1697 groups were identi- 
fied, and their magnitude gaps were compared with galaxy groups 
selected from the Millennium simulation at redshift z ~ 0.041. 
The magnitude gap distributions from both SDSS-DR4 as well as 
the Millennium simulation are shown in Fig.QT] Panels a and b of 
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Table 2. Comparison between the observed and simulated fraction of 
groups with magnitude gaps Amu and Amy, estimated from the his- 
tograms presented in Fig. 1 1 1 1 




3.5 3.0 2.5 2.0 1.5 
Look-back time (Gyr) 



1.0 0.5 0.0 



Figure 10. The fractions of F12 (black line) and F14 (thick red line) groups, 
identified at redshift z = 0, that survive as fossils, as a function of look-back 
time. 




Figure 11. The il-band magnitude gap distribution for haloes from the Mil- 
lennium semi-analytic models of Bo wer et al 1 fcOOfj) (red histograms) su- 
pe rposed on the r-b and data from Sample II of SDSS-DR4 group catalogue 
of lYang et alj 1 2007) (black histograms), (a) This shows the magnitude gap 
Amu between the the first and second most luminous galaxies, compared 
with galaxies from the SDSS-DR4 catalogue of groups computed within 
group radius, (b) This is the same as in (a), but for the magnitude gap Ami4 
between the first and the fourth most luminous galaxies. The 1697 SDSS- 
DR4 groups are within the mass range log(Af (i?igo)/^ — 1 Mq) ^ 13.25, 
and redshift range 0.01 ^ z < 0.1. Those from the Millennium simula- 
tions consist of 14612 X-ray groups selected at redshift 2 ~ 0.04 in the 
same mass range. 



Fig. [TT]refer to the magnitude gap distribution Amu and Amu 
respectively. Results show that the estimated magnitude gaps are 
in fair agreement with observation. The fraction of groups with 
Amy 55 2.0 is more or less the same while the fraction of those 
galaxy groups with Amu ^ 2.5 is different by ~ 1% (see Ta- 
ble. [2}- This shows that the incompleteness resulting from limited 
resolution does not affect our statistics. 



Selection criterion 


SDSS (DR4) 


Millennium simulation 




(Yang et al. 2007) 


(Bower et al. 2006) 


Amu >2.0 


2.0% ± 0.4 


2.1% ±0.2 


Ami4 >2.5 


6.2% ± 0.6 


5.1% ± 0.2 



6 DISCUSSION AND CONCLUSIONS 

In this work, we analysed the evolution of the magnitude gap (the 
difference in magnitude of the brightest and the nth brightest galax- 
ies) in galaxy groups. Using the Millennium dark matter simu- 
lations and associated semi-analytical galaxy catalogues and gas 
simulations, we investigated how the magnitude gap statistics are 
related to the history of mass assembly of the group, assessing 
whether its use as an age indicator is justified. 

A catalogue of galaxy groups, compiled from the Millennium 
dark matter simulations, was cross-correlated with catalogues re- 
sulting from hot gas simulations, and from semi-analytic galaxy 
evolution models based on these simulations. This resulted in a 
list of groups, with various properties of the associated dark mat- 
ter haloes and galaxies, at 21 time steps, over the redshift range 
z ~ 1.0 to z = 0. The simulated X-ray emitting hot IGM properties 
were known for these haloes only for z~0, and these were used to 
define a sample of X-ray emitting groups. This is necessary since 
our objective was to examine the evolution of fossil groups, which 
are observationally defined in terms of both optical and X-ray pa- 
rameters. 

We compared the estimated magnitude gaps in these galaxy 
group s fro m two different se mi-analytic models of iBower et al.l 
i feOOfj) and lCroton etall < l2006h . based on the Millennium dark mat- 
ter simulations, and found that the model of Bowe r et al. I d2006h 
better matches the observed present-day distribution of the differ- 
ence in magnitude between the brightest galaxy in each group, and 
the second and t hird brightes t galax ies Amu and Amu. We de- 
cided to use the lBower et al. (2006) catalogue for the rest of this 
study. 

We examined the evolution with time of fossil galaxy groups, 
conventionally defined as those with an .R-band difference in mag- 
nitude between the two brightest galaxies Amu ^ 2 (within 
0.5 -R200 of the group centre). We explored the nature of the groups 
that would be selected if the radius of the group were extended to 
R200, and the definition of the magnitude gap in terms of Amu 
were varied. Our major conclusions from the analyses can be sum- 
marised as follows: 

(i) The parameter Am^ defined for a galaxy system as the mag- 
nitude gap between the first and i th brightest galaxies (estimated 
within a radius of O.5R200 or R200) can be shown to be linked to 
the halo mass assembly of the system ai.o (Fig.[7J, such that galaxy 
systems with larger magnitude gaps Am^ are more likely to be 
early-formed than those with smaller magnitude gaps. 

(ii) Fig.|6]shows that, contrary to expectation, irrespective of the 
redshift at which fossil groups are identified according to the usual 
criteria, after ~4 Gyr, more than ~90% of them become non-fossils 
according to the magnitude gap criterion. Over the span of 7.7 Gyr, 
which is the time interval between z =0-1, very few groups retain 
a two-magnitude gap between the two brightest galaxies. This pro- 
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vides clear evidence that the fossil phase is a temporary phase in the 
life of fossil groups (also see lvon Benda-Beckmann et al.l 2008). 

(iii) In a given galaxy group, the merging of the i th brightest 
(or a brighter) galaxy, with the brightest galaxy in the group (often 
the central galaxy if there is one), results in an increase of Amu. 
However, one of the main reasons for the fossil phase to be a tran- 
sient one is that such a magnitude gap could be filled by the infall 
of equally massive galaxies into the core of the group, which would 
lead to a decrease in Ami,. Therefore, groups with smaller mag- 
nitude gaps are not necessarily late-formed systems. Many groups 
spend a part of their life in such a fossil phase, though an over- 
whelming majority of them would not fulfil the criteria of the "fos- 
sil" label at all epochs. 

(iv) For our sample of X-ray bright groups, the optical criterion 
Ami4 ^ 2.5 in the 7?-band is more efficient in identifying early- 
formed groups than the condition Amu ^ 2.0 (for the same filter), 
and is shown to identify at least 50% more early-formed groups. 
Furthermore, for the groups selected by the latter criterion, the fos- 
sil phase in general is seen on average to last ~ 1.0 Gyr more than 
their counterparts selected using the conventional criterion. 

(v) Groups selected according to Amy 2.5 at z — corre- 
spond to ~ 75% of those identified using the Amu 2.0 crite- 
rion. On comparing different panels in Fig. [8] one finds that early- 
formed groups identified from their large magnitude gaps (either 
Amu 2.0 or Amy ^ 2.5) represent a small fraction (18% 
for F14 and 8% for F12) of the overall population of early-formed 
systems. This is especially noticeable in the high-mass regime. 

(vi) Finally, Fig. [9] shows that in comparison to conventional 
fossils (i.e. Fia groups), the F14 groups identified based on 
Amy 2.5, predominantly correspond to systems with halo 
masses M{R 20 o) s= W 14 h- 1 M (D . 

This makes the criterion Amu 2.0 marginally more efficient 
than Amy ^ 2.5 in identifying massive early-formed systems. 



These results depend to some extent on the employed semi- 
analytic model in our current analysis, and the statistics might 
change if one uses different semi-analytical model of galaxy for- 
mation. 

Physical prescriptions such as galaxy merging, supernova and 
AGN feedback used in such models are somewhat different from 
one another. Furthermore, superfluous mergers may result from the 
algorithm used for the identification of haloes in the Millennium 
DM simulation, and this may affect the merger rates calculated 
from various stud ies, including ours, using these catalogues (e.g. 
Genel et al. 2009). Though merging is the most important process 
affecting galaxies in groups, there are other physical processes such 
as ram pressure stripping, interactions and harassment, group tidal 
field, and gas loss, that are not fully characterised by current semi- 
analytic models. 

This is partially due to the limited spatial resolution of the 
Millennium simulation. The new release of the current simulation, 
i.e. Millennium-II Simulation might help to address some of the 
above issues. The latter has 5 ti mes better spatial resolutio n and 
125 times better mass resolution dBovlan-Kolchin et al . 2009). Fu- 
ture semi-analytic models based upon high resolution simulations, 
incorporating such effects would be worth employing in a similar 
investigation to find better observational indicators of the ages of 
galaxy groups and clusters. 
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